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Abstract: We make a comprehensive investigation of the observational effect of the in- 
flation consistency relation. We focus on the general single-field inflation model with the 
consistency relation r = —8c s nt, and investigate the observational constraints of sound 
speed c s by using the Seven- Year WMAP data, the BICEP tensor power spectrum data, 
and the constraints on /nl" 1 ' an d /nl* 1 from the Five- Year WMAP observations. We find 
that the constraints on the tensor-to-scalar ratio r is much tighter if c s is small, since a 
large tilt nt is strongly constrained by the observations. We obtain r < 0.37, 0.27 and 
0.09 (dn s /dlnk = 0) for c s =l, 0.1 and 0.01 models at 95.4% confidence level. When tak- 
ing smaller values of c s , the positive correlation between r and n s also leads to slightly 
tighter constraint on the upper bound of n s , while the running of scalar spectral index 
dn s /dlnk is generally unaffected. For the sound speed c s , it is not well constrained if 
only the CMB power spectrum data is used, while the constraints are obtainable by taking 
/nl 111 ano - /ml* 1 ' P r i° rs i n to account. With the constraining data of f^ lL and /§L h ', we 
find that, c s < 0.01 region is excluded at 99.7% CL, and the c s = 1 case (the single-field 
slow-roll inflation) is slightly disfavored at 68.3% CL. In addition, the inclusion of /nl" 1 
and /nl 11 ' m ^ j0 the analysis can improve the constraints on r and n s . We further discuss 
the implications of our constraints on the test of inflation models. 
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1. Introduction 

An important task of modern cosmology is to understand the expansion history of the 
Universe. The standard hot big-bang model is successful in explaining various observa- 
tions, including Hubble expansion, Big-bang Nucleosynthesis and microwave background 
radiation Q, yet still suffers from the flatness, horizon and monopole problems, etc. The 
inflation model, in which the vacuum energy drives the Universe exponentially expanding 
in the very early Universe f2], was proposed under such concerns. Besides the successful 
explanation of the above problem, inflationary cosmology can provide a viable mechanism 
for the origin of the cosmic structures. 

There have been numerous inflation models proposed in the last several decades. In 
the face of so many competing candidates, it is necessary to find an effective way to figure 
out which one is realistic, or at least, which one is most favored by the cosmological 
observations. Especially, it is important to confirm or rule out the canonical single-field 
slow-roll (SFSR) inflation model. 

It has been proved that the SFSR inflation can generate observable primordial scalar 
and tensor perturbations, which encode themselves in the cosmic microwave background 
(CMB) anisotropics. Thus, it is possible to test inflationary models from the current CMB 
observations, e.g., the Wilkinson Microwave Anisotropy Probe ( WMAP) satellite Q, QUaD 
experiment [Q], BICEP experiment |J and other probes ||. 
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There have been a number of investigations made on testing the inflation models 
from the current and future CMB observations , mainly on the issues of constraining the 
SFSR model with the scalar spectral index n s , running of the spectral index dn s /d\a k, and 
tensor-to-scalar ratio r as free parameters. Besides the determination of the parameters in 
inflation models, it has also been proposed j8j that the consistency relations, which features 
various types of inflation models, can be used as a test to classify and distinguish different 
models of inflation. The possibility of the observational test of the consistency relations 
has been discussed in || in detail. 

In this paper we make further investigations on the observational effect of the consis- 
tency relation r = —8c s nt. We focus on the general single-field inflation model, and discuss 
the current constraints on the sound speed c s from the CMB data, including the Seven- 
Year WMAP (WMAP7) power spectrum data [§, the BICEP data §], and the constraints 
on f^ iL and /£ r L th - from the Five- Year WMAP (WMAP5) observations 0, We then 
discuss the results of the constraints on n s , r, dn s /dlnk parameters when c s 7^ 1. 

This paper is organized as follows. In Sec. 2, we introduce the inflationary consistency 
relation in the general single-filed inflation model. In Sec. 3, we briefly introduce the 
CMB data and data analysis methodology used in this paper. The results of constraints 
on cosmological parameters are presented in Sec. 4 and Sec. 5. We summarize our results 
in Sec. 6. 

2. Single- field Inflation Model 

Let's start with the general single-field inflation model 



S = I d A x^—g 



M l 



(2.1) 



where M p = 1 / V 8nG is the reduced Plack mass, R is the Ricci scalar, g is the determinant 
of the metric, and X = —-^g^ v d^<j)d v (f). P(X,(f>) is an arbitrary function of X and <ft. This 
action is the most general Lorentz invariant action for inflaton eft minimally coupled to 
Einstein gravity. The primordial scalar power spectrum of curvature perturbation is [12] 

A « = tp£ • (2 - 2) 



where 



is the slow-roll parameter, and 



(2-3) 



Px (2.4) 



P x + 2XP XX 

is the speed of sound. The spectral index of scalar curvature perturbation power spectra 
becomes 

dlnA 2 R . , 

n '- 1 = liMF = - 2t -"-'- < 2 ' 5) 
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where 



are another two slow-roll parameters. Due to the dynamics of inflation, the spectral index 
n s can be scale-dependent as well. Its scale-dependence is measured by the running of the 
spectral index dn s /d\n k. The primordial power spectrum of scalar curvature perturbation 
then takes the form, 

(t \ n s (ko)-l+^dn 3 /dlnk 

where ko is the pivot scale. The primordial power spectrum of gravitational waves pertur- 
bation generated during inflation only depends on the Hubble parameter during inflation 

2 _ H*/M$ 



with tilt 



din AS, , 

"* = h^f = - 2f - (2 - 9 > 

The tensor-to-scalar ratio is defined as 

r = A^/A| = 16c s e, (2.10) 

and then by combing with Eq. ( |2.9[ ), we obtain the consistency relation 

r = -8c s n t . (2.11) 

Here we should note that since the acceleration of scale factor takes the form a = H 2 a(l—e), 
inflation only happens if e < 1, therefore from Eqs. (|2.9|) and ( 2.1 1| ) , we know the valid 
ranges of values for n t and r as 

- 2 < n t < 0, and r < 16c s . (2.12) 

If c s = 1, the general single- field inflation model reduces to the single- field slow-roll 
inflation. But if c s <C 1, the non-trivial sound speed of inflation can generate non-Gaussian 
modes of perturbation, which results in a large non-local form of bispectrum [13|. Although 
the non-local form of bispectrum has not been well classified, the two most general types, 
equilateral type with shape size /nl"' > an d orthogonal type measured by f^ h ', have been 
widely discussed in literatures || [ll], [lj|]. In |[4|], the observational constraint on the c s 
from the bispectrum has been discussed and the requirement from the stability of the field 
theory (c 2 s > 0) implies /£ r L th - < -0.054/^ uiL - 
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3. Data Analysis Methodology 



In the following data analysis, we will combine WMAP7 power spectrum ||, with BICEP 
tensor power spectrum data || and bispectrum constraints on f^ lL and f^-^ 1 ' [|H]> > t° 
constrain inflation consistency relation. The WMAP TT power spectrum at 2 < / < 1200 
is powerful to constrain the cosmological parameters, e.g. n s and dn s /dlnk. We also 
use the WMAP TE/EE data at 2 < / < 800, and the BB data mainly on large scales 
2 < I < 23. To be consistent with the WMAP results || , we choose our pivot scale to be 
k = 0.002 Mpc" 1 . 

We also use the BICEP tensor power spectrum data which mainly covers the region 
21 < / < 335. Following the pipelines of ||, we construct the expected bandpowers for the 
inflation models, and use the lognormal approximation to calculate the \ 2 function, 

X 2 (P) = [z BB - Z( P ) BB ] T [D^(p)]- 1 [z BB - Z( P ) BB ] , (3.1) 

where p is the model parameters, and Z BB and Z(p) BB are the observational and the- 
oretical bandpowers. T) BB (p) is the covariance matrix which is dependent on the model 
parameters. The likelihood function then takes the form 

C oc 1 p -x 2 (p)/2 (3 . 2) 

v / det[D^(p)] 

In addition, the observational results of /jS ' ' and f^ h ' can also constrain the value 
of c s . We use /m, an d /nl^ priors obtained from the WMAP5 observations [10, 11], and 
construct the x 2 function as 

X 2 (P) = w(p)wMApC~Vp)wMAP, (3.3) 



where C is the covariance matrix given in [11|, and u(p)wmap is the difference between 



the observed and model values of f^ 31 ' and /^L h t' ' 



r(n > | (/nl Ui1 '(p)) - (/nl U1 ')wmap . , , 

W(P)WMAP - 1 ,i 0lth .,s X /?orth^ I' y 6A > 



(/nl '(p)) ~ (/nl Owmap 



The WMAP5 data yields to Q 1 



/X iL = 155±140, /§ r L th - = -149±110. (3.5) 

We will determine the best-fit parameters and the 68.3% and 95.4% confidence level 
(CL) ranges by using the Monte Carlo Markov chain (MCMC) technique. The whole set 
of our free parameters is 

P = {tt b h 2 ,£l c h 2 ,9,T,n s ,dn s /dlnk,r,c s ,A s ,Asz} 2 - (3.6) 



1 The covariance matrix C is dependent on the data jo]. Since the WMAP7 covariance matrix C has 
not yet been published, we will adopt the WMAP5 covariance matrix C in the following discussion. 

2 6 is the ratio of the sound horizon to the angular diameter distance; r is the the reionization optical 
depth; A a is the primordial superhorizon power in the curvature perturbation on the pivot scale ko = 
0.002 Mpc" 1 ; Asz is an SZ template normalization. 
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Figure 1: Power spectra for inflation models with different c s . Models with c s = l, 0.1 and 0.01 are 
plotted in blue solid, red dashed and green dotted lines, respectively. In all figures we fix r = 0.15. 
The primordial tensor power spectrum and the BB, TT, EE, TE power spectra are plotted. The 
WMAP data Q are plotted in black points. The c s =0.1 model with n t — —0.1875 leads to slightly 
larger values of power spectrum at the large scale (not very evident), while the c s =0.01 model leads 
to significantly larger Pr(k)/Cis in small-fc/low-Z region. 



We modify the publicly available CAMB |]| and COSMOMC packages fL6| to include 
models with c s as a free parameter, and generate O(10 5 ) samples for each set of results 
presented in this paper. 



4. Cosmological Constraints of Fixed c s Models 

In the following sections we will discuss the cosmological interpretations of the consistency 
relation. In this section we focus on three models with c s fixed as 1, 0.1 and 0.01. The free 
c s model will be discussed in the next section. 

4.1 Effects of c s on the Power Spectrum 

We firstly clarify how the different values of c s affect the shape of the angular power spectra 
of CMB. 

In the upper-left panel of Fig. [l] we plot the primordial tensor power spectrum for 
c s =l, 0.1 and 0.01 respectively. We see that c s has significant influence on the tilt of the 
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power spectrum through the consistency relation nt = — r/(8e s ). In particular, at the large 
scale (small k), the c s =0.01 and c s =0.1 models have much larger values of Pt(k) than that 
in the model with c s = 1. 

This effect is also visible in BB, TT, EE and TE power spectra, which is shown in 
the upper-right and lower panels of Fig. |l| (with lensing). In all figures we take r = 0.15 
and fix other parameters at their WMAP7 best-fit values. It is shown that the amplitude 
of the power spectrum for c<j=0.1 is slightly larger than the c s =l case (not evident), while 
for the c s =0.01 case the low-/ C;s are significantly larger. The panels indicate that the set 
of parameters r = 0.15, c s = 0.01 is inconsistent with the WMAP data. Thus, we expect 
a tight constraint on r when c s is small. 

4.2 Results of Fitting 

Our results of fitting for different models with fixed c s are shown in Table |l[ 



Table 1: Results of fitting with fixed c s . 



Data & Model 




r(95.4% CL) 


n t (95.4% CL) 


n s 


dris/dlnk 




Cs 


= 1 


< 


0.37 


> 


-0.05 


u - ao '-o.oio 




WMAP 


Cs 


= 0.1 


< 


0.26 


> 


-0.33 


972+0-016 
u - yl -0.014 






Cs 


= 0.01 


< 


0.09 


> 


-1.16 


966+ - 017 






Cs 


= 1 


< 


0.32 


> 


-0.04 


966+ - 024 




WMAP+BICEP 


Cs 


= 0.1 


< 


0.26 


> 


-0.33 


971+0-019 
u,t " -0.014 






Cs 


= 0.01 


< 


0.09 


> 


-1.16 


n 969+ 014 
u.»uy_o .013 






Cs 


= 1 


< 


0.36 


> 


-0.05 


^'^-o.oss 


-0 023+ ' 018 
U.UZO_ Q 024 


WMAP+BICEP 


Cs 


= 0.1 


< 


0.35 


> 


-0.44 




-0 022+ ' 021 
u.uzz_ 034 


(+dn s /d\nk) 


Cs 


= 0.01 


< 


0.10 


> 


-1.21 


1 n05+ 068 
i.uuo_ 013 


-0 019+ ' 005 



Here we list the results for c s =l, 0.1 and 0.01. In the 4-9 rows, the constraints on 
r, nt, n s and dn s /dlnk by using the WM^P+BICEP data are listed, divided into the 
dn s /d\nk = and dn s /d\nk ^ cases. For comparison, in the 1-3 rows we also list 
the results obtained by using the WMAP data alone with dn s /dlnk = 0. We discuss the 
results of constraints in Sec. 4.2.1 and 4.2.2 in detail. 

4.2.1 The dn s /d\nk = case 

In this subsection we briefly discuss the fitting results of models without including dn s /d In k 
as a free parameter. 

Let us first see the constraints on the tensor-to-scalar ratio r which determines the 
amplitude of the tensor power spectrum. The fitting results of r are listed in the second 
column of Table [l] and the likelihood functions are plotted in the left panel of Fig. ^. Using 
the WMAP data alone, we find a 95.4% CL constraint r < 0.37, which is well consistent 
with the result obtained by the WMAP 7-yr data (r < 0.36). As expected, we find the 
constraint on r becomes much tighter if the value of c s becomes smaller (see the left panel 
of Fig. |). Using the WMAP+BICVP data, we find r < 0.32, 0.26 and 0.09 (95.4% CL) 
for c s = 1, 0.1 and 0.01, respectively. This result is quite reasonable. According to the 
consistency relation n t = — r/(8c s ), if c s is small, a large r would lead to a large nt, leading 
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Figure 2: Fitting results for the fixed c s models with dn s /d\nk = 0. In the left and middle panels 
we plot the likelihoods of r and n s for the c s = 1, 0.1, 0.01 models in blue, red, green lines, and 
the results obtained by using the WMAP and WMAP+BIGEP are shown in dotted and solid lines, 
respectively. In the right panel we also plot the marginalized 68.3% and 95.4% CL contours in the 
r — n s plane. We see that smaller values of c s lead to tighter constraints on r, and thus tighter upper 
bound constraints on n s , due to their correlation. The inclusion of BICEP data slightly improves 
the constraints on r and n s for the c s = 1 case. 



to a large tensor mode on superhorizon scale, which is strongly constrained by low-Z CMB 



data. The necessary condition for inflation Eq. (2.12) is automatically satisfied by the 
constraint. The r = 0.15, c s = 0.01 case shown in the Fig. |l] is excluded. 

Except for r, another interesting issue is the fitting results of the scalar spectral index 
n s . Using the WMAP+BICEP data, we find n s = 0.966+J];^, 0.971±g;g^ and 0.969jlg$| 
(68.3% CL) for c s = 1, 0.1 and 0.01. All the results are consistent with Harrison-Zeldovich 
spectrum (n s = 1) at 68.3% CL. The likelihoods of n s in different cases are plotted in 
the middle panel of Fig. ||. Similar to r, we find that the upper bound constraint on n s 
also becomes tighter when c s is smaller. This effect is caused by the positive correlation 
between r and n s . In the right panel of Fig. ^, we plot the marginalized contours in the 
r — n s plane, which shows that the smaller c s is, the tighter constraints on n s and r. 

The results of constraints from BICEP data are as follows. For the c s =l model, we get 
r < 0.37(0.32) (95.4% CL) by using the WMAP{ WMAP+BICEP) data. The inclusion of 
BICEP data slightly improves the constraint by ~ 14%. However, for the c a = 0.1 and 0.01 
cases, since the constraints on r mainly come from constraints on rit by the low-Z WMAP 
data, the inclusion of BICEP data does not lead to significant improvement in the result 3 . 
The BICEP data can only affect n s through its correlation between r, so it does not have 
significant effect on the results of n s . 

4.2.2 The dn s /d\nk / case 

We now discuss the results of constraints with running of spectral index dn s /dlnk as a 
free parameter. 



3 The Two- Year BICEP data, which maps only ~ 2% of the sky, measures limited modes of perturbations, 
therefore does not contribute too much on the total constraining. 
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Figure 3: Fitting results for the fixed c s models with dn s /din k ^ 0. The marginalized likelihoods 
of dn s /dlnk and n s are shown in the left and middle panels. In the right panel we plot the 
marginalized 68.3% and 95.4% CL contours in the dn a /d\nk — n s plane. The influence of c s on 
dn s /dlnk is not very significant (left panel). The inclusion of dn s /dh\k significantly amplifies the 
allowed range of n s and shifts its central value from below 1 to above 1 (middle panel), since they 
are strongly anti-correlated (right panel; see also Eq. (2.7)). 



The fitting results are shown in the last three rows of Table [T] and Fig. |3[ For the three 
models we find similar constraints on dn s /dlnk, thus c s does not have too much influence 
on dn s /dln k. 

The most striking effect of the inclusion of dn s /dlnk is the significant amplification 
of the allowed region of n s . In the middle panel of Fig. [3|, we plot the likelihood functions 
of n s for dn s /dlnk = (solid) and dn s /d\nk ^ (dotted). We see that the width of n s 
likelihood is increased by nearly a factor of two. The inclusion of dn s /dlnk also changes 
the central values of n s from below 1 to above 1. These phenomena are caused by the 
strong anti-correlation between dn s /dlnk and n s (see the right panel of Fig. [3]). 

The inclusion of dn s /dlnk slightly releases the upper bound constraints of r. For 
the c s = 1, 0.1, 0.01 models, the 95.4% CL upper bounds on r are 0.32, 0.26, 0.09 for 
dn s /dlnk = and 0.36, 0.35, 0.10 for dn s /d\nk ^ 0. 



5. Cosmological Constraints of Free c s Models 

In this section we consider the more general case, i.e., treating c s as a free parameter. 4 We 
not only use the WMAP and BICEP power spectrum data, but also take f^ 11 ' and f^ 1 ' 
into consideration in order to constrain c s . A summary of the fitting results, including c s , 
r, m, n s and dn s /d\nk are given in Table ||. Notice that the word 'UCON' represents for 
'unconstrained'. 

5.1 Results of fitting Without /nl Prior 

In this subsection we discuss the results obtained by WMAP+BICEP data, without adding 
4 We sample c B in the range of (0.001,1). 
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Figure 4: Marginalized likelihoods of c s , r, dn s /dlnk and n s for the free c s case, obtained by 
using the WMAP and WMAP+BICEP data. Upper-left: The current CMB power spectrum data 
alone can not constrain c s . Upper-right: The inclusion of dn s /d\nk slightly widen the distribution 
of r. Lower- left: The inclusion of BICEP data does not affect the distribution of dn s /dlnk too 
much. Lower-right: The inclusion of dn s /d\nk greatly boardens the width of n s likelihood, and 
shifts its central value to be greater than unity. 



Table 2: Results of fitting with c 8 as a free parameter. 



Data & Model 


Cs 


r(95.4% CL) 


ret (95.4% CL) 




dn s /dlnk 


WMAP 


+dn a /d In k 


UCON 
UCON 


< 0.36 

< 0.46 


> -1.33 

> -1.44 


966+ 024 

-, n97 +0.065 
i - uz ' -0.051 


-0 026+ 023 


WMAP+BICEP 


+dn a /d In k 


UCON 
UCON 


< 0.32 

< 0.41 


> -1.30 

> -1.40 


n q fi7 +0.025 

u.yof_ 012 

I nin+0.063 
i - ui3 -0.03S 


-0 026+ - 020 

U.UZD_ () Q27 


WMj4P+BICEP+/nl 


+dn s /dln k 


019+ 0012 

016+ 017 
u.uxo_ 003 


< 0.21 

< 0.29 


> -0.91 

> -1.15 


Q 973+0.011 

1 016+ 064 
i.uio_ 045 


-0 024+ - 022 
u - uz ^-0.033 



/ N l" an d /§l P r iors into the analysis. The fitting results of parameters are shown in the 
1-4 rows of Table |2|, and the likelihoods of c s , r, n s and dn s /d\nk are plotted in Fig. ^. 

Let us first have a look at the constraint on c s , and its likelihood is plotted in the 
upper-left panel of Fig. |j. The current CMB power spectrum data is not able to constrain 
on c s , and the likelihood function shows that it can take any possible value given the 
current constraints. 
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Figure 5: Marginalized 68.3% and 95.4% CL contours in the r — n t (left) and c s — n t (right) 
planes which are obtained by using the WMAP+BICEP data. In the left panel, the nt = — § line 
is plotted in the black dashed line. We see that n t > —2 automatically satisfied, and the inclusion 
of dn s /dhik evidently amplifies the parameter space. 



Secondly, the likelihoods of r are shown in the upper-right panel of Fig. ^. We see 
the green and blue lines are close to each other, which means that the constraint on r 
with free c s is similar to the result for c s =l. 5 We get r < 0.32 (95.4% CL) for both 
c s = 1 and c s free models ( WMAP+BICEP, dn s /d\a k = 0). The inclusion of BICEP data 
slightly improves constraint of r from 0.36 to 0.32 {dn s /d In k = 0). We see the inclusion of 
dn s /dlnk as a free parameter boardens the width of distribution of r. At 95.4% CL, the 
constraint is widened from r < 0.32 to r < 0.41 ( WMAP+BICEP). 

Thirdly, in the lower panels of Fig. || we plot the marginalized likelihoods of dn s /din k 
(left) and n s (right). We find that the results are similar to the fixed c s models. It implies 
that the BICEP data almost does not affect the constraint of dn s /dlnk, and the inclusion 
of dn s /dlnk greatly amplifies the distribution of n s and shifts its central value to above 1. 

In addition, in Fig. [| we plot the marginalized 2D-contours in the r — nt (left) and 
c s ~ n t (right) planes. The slow-roll model with nt = —r/8 is plotted in the black dashed 
line in the left panel, and the dn s /dlnk = and dn s /dliak ^ cases are shown in orange 
dotted and black solid lines, respectively. By releasing c s as a free parameter, nt ranges 
to much smaller values (~-l to -1.5 at 95.4% CL), especially in the small c s region. We 
find nt > —2 is still satisfied by the constraints. In both panels we see the inclusion of 
dn s /dh\k boardens the ranges of the parameter space. 



5.2 Results of fitting With /nl Prior 

In this subsection let us take the f^ 11 ' and f^ h ' priors into account. The fitting results 
are shown in the last two rows of Table ^ and the likelihoods of c s , r, n s and dn s /dlnk 
are plotted in Fig. |(| 



3 A similar conclusion was obtained in 
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Figure 6: Marginalized likelihoods of c s , r, n s and dns/dlnk for the c s free case from 
WMAP+BICEP data. The inclusion of /nl prior evidently tightens the constraints on c s and 
r, and slightly improves the upper-bound constraint on n s . The running of the spectral index 
dn s /dhik remains unchanged. 




c s free, WMAP+BICEP 
— With/v,. 
— - Without /vl 



■a 

5 




c, free, WMAP+BICEP 

— With/m. 

— Without/^ 

I I 



0.95 



1 



1.05 
n. 



1.2 



Figure 7: Marginalized 68.3% and 95.4% CL contours in the r — n s (left) and n s — dn s /dlnk 
(right) planes, obtained by using the WMAP+BICEP data. In all figures we let dn s /dlnk and c s 
as free parameters. The constraints with and without /nl cases are shown in dark cyan solid and 
orange dotted colors respectively. 



The most striking effect is that the inclusion of f^ lL and f^^ 1 ' significantly improves 
the constraint on c s . At 68.3% CL, we obtain 0.013 < c s < 0.031 and 0.013 < c s < 0.033 
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for the cases of without and with dn s /d\nk as a free parameter. We find the large c s 
region, including the SFSR inflation with c s = 1, is slightly disfavored at around 68.3% 
CL, while c s < 0.01 is excluded at 99.7% CL. Thus, the bispectrum data is much more 
powerful than the power spectrum data for constraining c s . 

By narrowing the allowed range of c s , the addition of f™f* 11 ' an d /nl* 1 a ^ so nas 
teresting effect on constraining the other parameters. The likelihood of r is shown in the 
upper-right panel of Fig. ^. We see that, once /nl"' an< ^ fm^' P r iors are considered, the 
constraint becomes much tighter. This can be also seen through the contours in the left 
panel of Fig. [5]. Again, due to the correlation between r and n s , where /nl" 1 anc ^ /nl^ 
priors are included a slightly tighter constraint on the upper bound of n s is also obtained 
(see the lower-left panel of Fig. ^). But the effect of f^jf' an< ^ /nl^ P r i° rs on dn s /dlnk 
is negligible (see the lower-right panels of Fig. |6|) . 

Finally, the correlations between r, n s and dn s /dlnk are shown in Fig. [7|. The left 
panel shows the r — n s contours and the right panel shows the n s — dn s /dlnk contours. In 
all figures we set c s , dn s /dlnk as free parameters and use both WMAP and BICEP data. 
The cases of without and with /nl as a free parameter are shown in orange and dark cyan 
colors. One can see the strong correlation between r and n s , which suggests that if the 
distribution of r is tightened, n s distribution is also constrained. However, the distribution 
of running spectral index dn s /dlnk, is not much affected by this correlation, because the 
change of n s is much smaller comparing with r. 



6. Conclusion 



In this paper we make a detailed investigation of the cosmological interpretation of the 
consistency relation from CMB data. We focus on the general single-field inflation model 
in which the spectral index rit of tensor perturbation power spectrum is related to the 
tensor-to-scalar ratio r by nt = — r/(8c s ), and further investigate the effect of the sound 
speed c s . The datasets used in this paper include the WMAP power spectrum data, the 
BICEP B-mode polarization data, and f^^' and /§L h ' priors obtained from the WMAP5 
bispectrum data. 

We discuss three models with fixed c s =l, 0.1 and 0.01. We find that when c s is small, 
the tilt of the tensor power spectrum nt becomes very large if r is not too small, and then a 
tight constraint on r is obtained for c 8 <l. Using the W^M^P+BICEP data, we obtain the 
95.4% CL constraints of r <0.37, 0.26, 0.09 for the c s =l, 0.1, 0.01 cases (dn s /dhxk = 0). 
Due to the positive correlation between r and n s , smaller values of c s lead to slightly tighter 
constraint on the upper bound of n s . The effect of c s on the running of scalar spectral index 
dn s /d\nk is not obvious. The inclusion of dn s /dlnk significantly alters the constraints of 
n s , and slightly amplifies the upper bound constraint of r. 

For more general cases in which c s is taken as a free parameter, we find that c s un- 
constrained if we only use the current CMB power spectrum data in the analysis, and 
the marginalized distribution of r, n s and dn s /dlnk are all similar to the c s = 1 case. 
However, after taking f^ lh and f^ 1 ' priors into consideration, we find the sound speed 
c s is effectively constrained. The c s < 0.01 region is ruled out, and the c s > 0.03 region 
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is disfavored at the 68.3% CL. Prom the constraints on c s , the inclusion of /nl leads to 
tighter constraint on the r and n s . In the dn s /d\nk = case, we find r < 0.21/0.32 (95.4% 
CL) with/without /nl prior {dn 8 /d\nk = 0), and the results for the dn s /dlnk / case is 
r < 0.29/0.41. The running of spectral index dn s /d\nk is almost unaffected. 

To summarize, we find that the consistency relation has significant effect in the con- 
straints on cosmological parameters r and n s when c s is small, while the parameter 
dn s /dlnk remains unaffected. Although the sound speed c s is unconstrained by the CMB 
power spectrum data, it can be effectively constrained by the CMB bispectrum data. Using 
the /nl 111 ano - /nl 11 ' P r iors obtained from the WMAP5 data, we find the SFSR model with 
c s =l is slightly disfavored at the 68.3% CL. Thus, we are expecting that, the on-going and 
upcoming CMB observations, such as Planck [18| and CMBPol fl~9[l , with much lower in- 
strumental noise and better foreground clean, will provide stronger constraints on inflation 
models. 
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